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C^ . Abstract 

We study in detail the exact thermodynamics of the one-dimensional standard 
and staggered spin-1 Ising models with a single- ion anisotropy term in the presence 
of a longitudinal magnetic field at low temperatures. The results are valid for the 
'^ \ ferromagnetic and anti-ferromagnetic (AF) models and for positive and negative 

G ■ values of the crystal field for T > 0. Although the excited states of the ferro and 

anti-ferro models are highly degenerate, we show that the temperature required 
for reaching the first excited state in the classical spin-1 ferro model gives a scale 
of temperature that permits fitting the z-component of the magnetization only by 
K»- ■ the contribution of two ground states of the model. This approximation is not 

true for the equivalent AF function due to the fact that the AF model is gapless 
along the lines separating the phases in its phase diagram at T = 0. We relate the 
•n ■ number of plateaus in the magnetization of each model to their respective phase 

t:!" , diagrams at T = 0. The specific heat per site of the AF model distinguishes, at 

^^ I low temperature, the transitions A f^ E and G f^ E as the external magnetic field 

is varied. The exact Helmholtz free energy of the classical spin-1 model is mapped 
onto the equivalent function of the ionic limit of the ID extended Hubbard model 
by proper transformations. 
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1 Introduction 

Exactly solvable models help us having insights on more complex systems. In one space 
dimension (1-D), the thermodynamics of the spin-| Ising model in the presence of a 
longitudinal magnetic field has been solved in the whole interval of temperature using the 
density matrix approach[Tl |21 E]- To do so, it is essential the commutative nature of all 
operators in its Hamiltonian. This model is usually called a classical system. 

Cold atoms have made possible simulations of spin models. Recently Simon et al. [1] 
simulated a one-dimensional spin-1/2 model in the presence of a magnetic field with 
longitudinal and transverse components by using a Mott insulator of spinless bosons in 
a tilted optical lattice. The device was used to study the phases of this classical model 
at low temperatures. These optical devices have permitted the experimental study of 
properties of chain models, including the spin models. 

The spin-1 Ising model with single-ion anisotropy term, the Blume[5]-Capel[6] model, 
in the presence of a longitudinal magnetic field, has also a classical nature. This fact 
permits us to apply the transfer matrix method [H [21 [3] to solve it. With this approach, 
Aydiner and Akyiiz [7| obtained the numerical solution of the ID spin-1 anti-ferromagnetic 
(AF) Ising model with a single-ion anisotropy term in the presence of a longitudinal 
magnetic field. They studied the magnetic and thermal behavior of the model at very low 
temperatures. In 2003 Chen et a/. [8] applied the classical Monte Carlo (CMC) technique to 
the numerical calculation of the phase diagram at T = of the spin-1 anti-ferromagnetic 
Ising model in the presence of a single-ion anisotropy term for positive crystal field. 
The CMC was applied also to the study of magnetization plateaus of this model at low 
temperature. 

In 2005 Mancini[9] used the mapping between the ID extended Hubbard model in 
the ionic limit (i.e., |t| <^ \U\, \V\, in which t is the hopping exchange, U is the on-site 
Coulomb interactions and V is the inter-site Coulomb interaction) and the spin-1 Ising 
model with a single-ion anisotropy term in the presence of a longitudinal magnetic field to 
write the exact Helmholtz free energy (HFE) of the latterj^] as a set of coupled equations. 
The hierarchy of the equations of motion is closed, and they can be solved numerically. 
In Ref. [To] they have extended the approach to include a biquadratic nearest-neighbour 
interaction term, numerically solving the coupled equations and presenting the behavior 
of some thermodynamic functions at low temperatures. We derived in Ref. [TT] the exact 
expression of the HFE of the ID spin-1 Ising model with the {SzY term in the absence of 
a magnetic field. This HFE has a simple analytic expression valid for the ferromagnetic 
and AF models in the whole range of temperature. The absence of a magnetic field in this 
solution prevents accessing some information about the system, e.g.: information about 
the magnetization and the susceptibility of the model. 

All the previously mentioned papers on the one- dimensional spin-1 Ising model seem 
to be unaware of the nice work by Krinsky and Furman[T2] in 1976, which presents the 
calculation of the exact HFE for the ID ferromagnetic spin-1 Ising model with a single-ion 
anisotropy term, a biquadratic nearest-neighbour interaction term and a nonsymmetric 
term, in the presence of a longitudinal magnetic field. More recently Litaiff et al. also 
applied the transfer-matrix technique to calculate the exact expression of the HFE of 
this model. The authors did not use this thermodynamic function to explore its low 



temperature behavior. 

Although the ID spin-S* Ising model with the single-ion anisotropy term in the pres- 
ence of a magnetic field is called a classical model, it has a quantum nature that mani- 
fests itself as plateaus in the z-component of the magnetization at very low temperature. 
The existence of those plateaus has been demonstrated in a number of ID mo dels p^ 
and experimentally measured in some materials [I5 | fT6 | fT7] described by ID spin-models. 
Magnetization plateaus in the the AF spin-1 Ising model with positive crystal field in the 
presence of a longitudinal magnetic field have been obtained numerically in Refs. [3 El HO] ■ 
The authors of Ref. [12] mentioned the discontinuous change of the magnetization of the 
ferromagnetic model, but their focus was on the exact renormalization group of the model 
and the discussion of its critical points. 

The study of the simple one-dimensional spin-1 Ising model helps to understand the 
origin of plateaus in the magnetization function. In the present communication we present 
the exact analytic expressions of the HFE's of the standard and the staggered versions of 
the ID spin-1 Ising model with the {Sf)"^ term in the presence of a longitudinal magnetic 
field, valid for T > in section HJ Our solutions apply equally well to the ferromagnetic 
and to the AF models of both versions, extending the validity of the solution derived in 
Ref. [T2] for the ferromagnetic model. We study the z-component of the magnetization 
and the entropy per site at very low temperatures for the ferro and AF versions of the 
classical spin-1 model in sections [3] and HI respectively. We relate the number of plateaus 
in the magnetization to the number of phases in the ground state of the model at T = 0. 
Finally in section [S] we present our conclusions. More detailed calculations are shown in 
Appendices \M and [Bl 

2 The Hamiltonians of the classical spin-1 models in 
the presence of a longitudinal magnetic field 

The Hamiltonian of the one- dimensional of the classical spin-1 Ising model with a single- 
ion anisotropy term, the Blume-Capel model, in the presence of a longitudinal magnetic 
field isia E] 

N 



a' = Y. [JS^sUi-h's^ + D'isty] 



i=l 

N 

y , Hj^j+i, (1) 



where S^ is the z- component of the spin-1 operator with norm: IJS'IP = 2. J is the 
exchange strength and it can have either a negative value (ferromagnetic model) or a 
positive value (AF model). The crystal field D can assume positive, negative or null 
values. The model satisfies spacial periodic boundary condition in a chain with N sites. 
The external magnetic field is oriented along the easy-axis z. In this paper we use natural 
units, e = m = h = 1. Comparing our Hamiltonian ([1]) to the Hamiltonian (25) in Ref. [13] 



we have J = J*, h' = jjbH* and D' = D*, in which J*, H* and D* are the parameters of 
the Hamiltonian of the single spin-1 Ising model in Ref.|13]. 

Mancini[9j showed that by the substitution Sf = rij — Ij (where Ij is the identity 
operator at the i-th. site, in which rij = nj_^ + rij^j,, rij^o- = cIo-Cj.o- with a G {t^i}, c]^ is 
the fermionic creation operator of an electron at the z-th site with the spin component a 
and Cj^o- is the corresponding destruction operator, the Hamiltonian ([1]) is mapped onto 
the ionic limit of the Hamiltonian of the ID extended Hubbard model in the presence of 
a chemical potential, that is 



Il' = IlHub + N{J + h' + D')l, (2) 

where millH] 

N 

Unub = ^ (f/rii^^n,,; + 2^11^ rii+i - fiUi) . (3) 

The last equality is valid for U = 2D', V = ^ and fj. = 2J + h' + D'. As a conse- 
quence of relation ([2]) we have, for these values of parameters, that WHub{U,V, fi; (3) = 
>Vi( J, h', D'; (3)-{J + h' + D'), where Wnub (>Vi) is the HFE of the extended Hubbard 
model in the ionic limit (the HFE of the spin-1 Ising model ([T])). One is reminded that 
(3 = -gf, where k is the Boltzmann's constant and T is the absolute temperature in kelvin. 

The Hamiltonian ([T]) in the presence of a staggered longitudinal magnetic field is 



2M 

HL, = E [Js^su^-h'i-iys^ + D'isn']. (4) 

The Hamiltonian continues to satisfy space periodic condition, but now we have an even 
number of sites in the chain: A^ = 2M. 

Applying the method presented in Ref.[2n] we calculated in Ref.[2I] the /3-expansion 
of the HFE of the one- dimensional spin-S" Ising model, with single-ion anisotropy term, 
in the presence of a longitudinal magnetic field, VVs^J, h', D'; (3), up to order /3^^, in the 
thermodynamic limit. The /3-expansion of the HFE of the staggered Hamiltonian (jlj), 
with D' = 0, Wg"'^{J, h', 0; /3), was derived in Ref.[22] up to the same order in /3 also for 
M — >■ oo. Unfortunately those expansions do not permit us to study the thermodynamics 
of the standard and the staggered ID spin-1 Ising model close to T = 0. In our web pagqj 
we provide the data files with the quantum (arbitrary spin-S*) and the classical HFE's of 
the normalized Hamiltonians ([T]) and (jlj) up to order /3^^ for both versions (standard and 
staggered) of the model. 

In the present work we apply the transfer matrix method [3] together with the /3- 
expansion of the function VVi( J, h', D'] /3) obtained from the results of Ref.|21j with 5 = 1 
to calculate the exact HFE of the Hamiltonian (Jlj) in the thermodynamic limit (A^ — )■ 



^ Oiir web page: http://www.proac.uff.br/mtt 



oo), valid for T > 0. The expression of Wl°'^{J,h',D';j3) is obtained by using a well 
known result in the literature [2^ , namely, Wf^'^i-lh', D'; /3) = Ws{-J,h', D'; /3), with 
5" = ^, 1, |, ■ • • . This equality is valid for T G [0, oo). 

In order to calculate the exact function yVi{J,h',D';P), valid for T > and in the 
limit of A^ — )■ oo, using the transfer matrix method, we rewrite Hamiltonian ([1]) as a 
symmetric operator in the i-th and {i + l)-th sites, that is, 

Hgl, = JS^SU, - hSt - hSl, + D{StY + D{Sl,Y. (5) 

By comparing eqs.([T]) and (E]), we verify that the /3-expansion of the HFE of this model, 
in terms of the parameters J, h and D in Hamiltonian (E]) may be obtained by applying 
the following change of variables: h' — t- 2h and D' — > 2D, in the expansion of Ref.pT] 
with 5* = 1. 

In appendix [A] we calculate the three roots of the third degree equation derived from 
the transfer matrix method for the classical spin-1 model with a single-ion anisotropy 
term in the presence of a longitudinal magnetic field. The root with the largest modulus 
gives the expression of the HFE of the model (11])/ ([5]) (see eq. flA.4|) ). 

In order to verify which root, eqs. flA.8a|) - flA.8c|l . corresponds to the eigenvalue Ai of 



matrix U, assumed to be the root with the largest modulus, we calculate the /3-expansion 
of the functions — 4ln[sj], i = 1,2,3, and compare each one with the expansion of the 
HFE of Ref.[2I] with 5* = 1 (having made the change of variables h' = 2h and D' = 2D in 
the /9-expansion of Ref. j21]). By direct comparison of the expansions, we obtain for finite 
value of /3, up to order /3^°, that the root si, eq. (IA.8a|l . is the eigenvalue Ai, that has the 



largest modulus among the eigenvalues of matrix U. Our analysis is valid for finite value 
of (3, which excludes the value T = 0. 

Our solution Si is equal to the eigenvalue of the cubic equation with the largest modulus 
calculated in Ref.[T2] for the ferromagnetic model (J < 0). The way we write the solutions 
(lA.9a|) - (]A.9c|) avoids the necessity of defining a cut in the complex plane for their 



calculation. Our results extended those derived in Ref . y^ to include the HFE of the AF 
spin-1 Ising model in the presence of a longitudinal magnetic field. 

The exact HFE of the standard one- dimensional spin-1 Ising model with the single-ion 
anisotropy term in the presence of a longitudinal magnetic field h is 



Wi{J,h,D-/3) = -^ In 



-Q cos(-) + — 



(6) 



in which P, 9 and Q are given respectively by eqs.( IA.7a|) . ( ]A.9a|) . ( ]A.9b|) . The HFE 



is an even function of h: Wi{J,h,D; (3) = Wi{J,—h,D; (3). One is reminded that 
Wf '^^(J, h', D'- /3) = >Vi(-J, h', D'- /3) and WnuhiU, V, /x; /3) = Wi(J, h' , D'; (3) - {J + h' + 
D'), with U = 2D', V = ^andfi = 2J+h' + D'. 

From the simple expression ([6]) one can calculate the thermodynamic functions of the 
standard and the staggered models, as well of the extended Hubbard model in the ionic 



limit in the absence of an external magnetic field, for any finite value of P, that does 
not include the temperature T = 0, that is, T G (0, oo). It is a well known fact that at 
T = 0, h = and D = 0, two eigenvalues of matrix U are degenerated. Our results are 
valid in the limit of T — )■ 0. Due to the presence of the longitudinal magnetic field in the 
HFE ([6]) we can derive from it the z-component of the magnetization and the magnetic 
susceptibility per site for the ferromagnetic (J < 0) and AF (J > 0) versions of the 
standard Hamiltonian (DO)/ dS]) as well as the staggered magnetization and the staggered 
magnetic susceptibility[22j of the Hamiltonian (^ for any value of the exchange strength 
J. For h = 0, the expression ([6]) of the HFE Wi( J, 0, D; (3) coincides with the exact result 
presented in Ref.[TT]. 

Our exact result permits us to study the plateaus of the standard and the staggered 
z component of the magnetization per site of the ferro and the AF models of the stan- 
dard and the staggered versions, respectively, of the one-dimensional spin-1 Ising model, 
with the single-ion anisotropy term, in the presence of a longitudinal magnetic field for 
temperatures close to T ~ 0. The authors of Ref. [12] mentioned the discontinuity of 
the magnetization and the magnetic quadrupolar moment of the ferromagnetic model at 
T = 0, not going any further on this point. Aydiner and Akyiiz [7] and Chen et al [8] stud- 
ied numerically these plateaus in the AF version of Hamiltonian (II])/(IS])- More recently, 
Mancini and Mancini[Tn] solved (also numerically) the self consistent equations that yield 
the exact HFE of the AF S = 1 Ising model with a biquadratic nearest-neighbour term 
in the presence of a longitudinal magnetic field. 

In the rest of this paper we let J = —1 for the ferromagnetic model and J = 1 for the 
AF model. The values of the parameters D, h and T are given in units of I Jl, that is: tt, , 

TjT and TjT^ respectively. 

We restrict our presentation mainly to the behavior of two thermodynamic 
functions associated to the Hamiltonians ([I])/® and (jl]) that helps us to 
understand the presence of plateaus in the classical spin-1 model: the z- 
component of the magnetization [21], Aiz{J-ih,D] (5) {M.z = —\ ^^)-i of the stan- 
dard model ([I])/ ([5]), and the staggered ^-component of the magnetization [22] . 

Mf''^{J,h,D;(3) (Mf""^ = i [{SI) - {SI)] = -i ^^^), of the staggered Hamiltonian 

(HD; and the entropy per site, S{J,h,D;(3) (s = P^ ^V where W = Wi(>Vf"^) for 

the standard (staggered) model. Those are studied at very low temperatures. The HFE's 
Wi and yVi "^ are even functions of h, therefore S{J, h, D; (3) and the magnetization func- 
tions (A^^( J, h, D] (3) and J^f'^^{J, h, D; /?)) are even and odd functions of h, respectively. 
For this reason we restrict ourselves to the case oi h > 0. 

In Fig{T] we present the phase diagram of the ferro (J = —1) and AF versions of 
Hamiltonian ([I])/([5]) at T = 0. In its caption we describe the meaning of the phases in each 
diagram. It is interesting to notice that the exchange coupling term J in the ferromagnetic 
Hamiltonian ([I])/([S]) favors the states with z component of the spin sf = =f1 and parallel 
neighboring spins. The same term in the AF version of this Hamiltonian also favors the 
states with s|^ = ^1 but for anti-parallel (Neel state) neighboring spins. In the ferro and 
AF models the Zeeman term favors states with sf = =f1 aligned with the external magnetic 
field h, whereas the single-ion anisotropy term shows two distinct behaviors: for D < 



states with sf = =f1 are favored, independently of their relative alignment; and for D > 
states with s^ = are favored, the spin being perpendicular to the external magnetic field 
applied. For D < all the terms in Hamiltonian ([I])/® force the neighboring spins to 
align to each other and consequently to the external magnetic field at T = 0. The ground 
state is a collective stable state under small temperature fluctuations, a point which will 
be made clear after we discuss the necessary temperature to excite the first excited state 
of the ferromagnetic chain. For D > the effect of the crystal field D competes with the 
exchange strength J and the external magnetic h. In the near future we will verify that 
in this case we need a smaller energy to break the alignment between neighboring spins 
and excite the first excited energy level of the ferromagnetic chain. 

For TjT > and only for Tjr > | the Fig{T^ has two distinct phases. The line between 

the phases is: tti = |7| ~ ^- In the region A of FigH^, at T = 0, the single-ion anisotropy 
term gives the main contribution to the ground state and the spins are perpendicular to 
the longitudinal magnetic field. Our phase diagram Fig{T^ coincides with Fig. 2 of Ref. [T2] 
with K = 0. 

Figjl}) presents the phase diagram of the standard AF model ( J = 1) at T = 0. For 
Iji < and A > there are two distinct phases that correspond to the competition 
between the exchange coupling term J and the Zeeman term. One of them is the phase G 
that is the Neel state when < tjt < 1; the exchange term gives the largest contribution 
to the ground state energy. For 7^7 > and tjt > the three terms in the AF Hamiltonian 
([I])/® favor distinct states at T = 0. Their competition is responsible for a richer phase 
diagram at T = for the AF model in this region of parameters. On the other hand, 
this competition among the terms of Hamiltonian ([I])/® makes the ground state of the 
one-dimensional system less stable under small fluctuations of the thermal energy, as we 
will verify in the discussion of the first excited state of the AF chain. The diagram of this 
model has two tricritical points: V{-rj-^ = 1 and D = 0) and Qirji = | and tj| = |) in 

FigJT]D. For < jji < I and tjt > 0, the AF model has three distinct phases at T = 0, 
namely: the G (Neel) phase, the E phase (in which half of the spins in the chain are 
perpendicular to the external magnetic field), and the B phase (in which all the spins are 
aligned with the magnetic field). We also have three phases for tjt > |, but in this region 
of the phase diagram at T = the Neel state is not one of the possible ground states 
of the system; there is an A phase instead, in which all spins are perpendicular to the 
longitudinal magnetic field, besides the phases B and E. The phase diagram FigJTjD agrees 
with Fig. 2b of Ref.[H] and shows qualitative agreement with Fig.2 of Ref. [7]. 

The transitions between the ground states of the standard AF model happen for the 



following values of magnetic field A > and the intervals of tjt : 

G — B : 1^ < and A = 1, (7a) 

G — E: 0<T^<| and — = 1-—, (7b) 

A^E: ^>1 and A = ^, (7c) 

h D 

E — B : 1^1 > and — = 1 + — . (7d) 

The phase diagram of the staggered ferromagnetic (J = —1) Hamiltonian (jl]) at T = 
is shown in FigJTJo whereas the diagram of the staggered AF (J = 1) model at T = is 
depicted in FigJT^. 

To the best of our knowledge, there is no detailed discussion in the current literature 
regarding the presence of plateaus in the z-component of the magnetization in ID ferro- 
magnetic spin models at very low temperatures, although its discontinuity is mentioned 
in Ref412j. From the exact expressions of the 2r-component of the magnetizations of the 
ferro and AF spin-1 classical models, valid in the thermodynamic limit, we verify that 
they have no discontinuity for T > 0. Rather, there is a range of values of the magnetic 
field for which a continuous transition from one magnetization plateau to another takes 
place. Our result ([6]) is equally valid for the ferro and AF models of the standard and 
staggered spin-1 Ising model with {SfY term in the presence of a longitudinal magnetic 
field T > 0. 

Our result (E]) of the HFE of the ferro (J = —1) and the AF (J = 1) Hamiltonians 
Q/dS]) are valid for -i^ > 0. A natural way to describe those models at finite tempera- 
ture is through their respective density matrix operator. From now on, we restrict our 
discussions to the condition tjt > 0. 

3 The ferromagnetic classical spin-1 model 

Let us consider the ferromagnetic model ([I])/ ([5]) with J = — 1. The ground states of the 
chain at T = in the phases A, B and C in FigH^ are 

|^o)a = |0)i®|0)2®---®|0)^, (8a) 

|^o)b = |l)i® |1)2®---® |l)jv, (8b) 

|^o)c = \-1)i®\-1)2<^---®\-1)n- (8c) 

where Sf\s)i = s\s)i, s = 0, ±1 and i = 1,2,- ■ ■ ,N. Their respective ground state energy 
are named Eq, Eq and E^ . The ground state (IHcj) is presented since at /i = its energy 
degenerates into the energy of the state |\E'o)a- 

The values of the ground state energies in units of \J\ are 



El 



0, (9a) 



E^ _ / 2h 2D\ 



^0 



^ -1 - ^ + ^ h (9b) 



1^1 V 1^1 \J\J 

where A^ is the number of sites in the chain. 

The density matrix operator of the ferromagnetic Hamiltonian ([I])/(jS]) is 

/ D h \ e-^^o e-^^" 

/^ -1' HT' HT' /^ = -^r- \^o)a a{^o\ + -^r- l^o)s s(^ol + 



e /3S? e-/3s- ^ 



f^l^o)cc(^o| + ^x5^ l$i'0(*r|---, (10) 



Zi ' "'^ ^^ "' Zi 



«=i 



where Zi is the partition function associated to the ferromagnetic Hamiltonian ([I])/©, 
Z\ = Tr[e~'^^], with (3 = -^- Ei is the first excited state of the chain and its degeneracy 
is equal to N for all three ferromagnetic phases. The first excited state is obtained 
from the ground state vectors fl5aj) - fl5c|) by flipping one of its spin-1 in the chain. In the 
thermodynamic limit, the first excited state is highly degenerate as well the higher excited 
states[25]. In principle, the contribution of all excited states of the chain should be taken 
into account for non null temperatures. We have no mathematically sound argument to 
affirm that, at low temperatures, the expansion (TTOj) should be cut after the contribution 
of the first excited state. The exact expressions of the thermodynamic functions can 
be derived from result (E]) of the HFE i^ of the HFE, and they show plateaus in the 
^-component of the magnetization of this model at low temperature. The comparison 
between the exact result and a proposed approximation expression of A^^ should be able 
to say how good is the latter. 

It is common sense that the the least value of temperature for which the contribution 
of the first excited state must be taken into account for the behavior of the thermodynamic 
functions is determined by the factor e~ ^t , But in the present model the excited states 
have degeneracy at least of order N , so all of them should contribute to the functions at 
non zero temperature. 

In the following discussion we determine in each phase of the diagram [1^ the value of 
temperature Tmax where 



^ -.,„.. _ e"^^ ^ 3.06 X 10"^ (11) 

where Eq is the ground state energy of the phase and Ei is its first excited state. Then 

kTmax El — Eq 

We want to verify if the temperature T^ax at each ferromagnetic phase is such that, for 
temperatures lower than Tmax, the dependence of the z-component of the magnetization 
on the external magnetic field has a "step-like" form. 

9 



We have three distinct first excited states of the ferro classical spin-1 model ([I])/®. 
They depend on the values of the parameters t4 and -Mr. 

There are two distinct first excited states in phase B (see FigH^): 
1.1) for A > and 1^ < -1 - A 

\^?)b = |l)i®|l)2®---®|-l)i®---®|l)jv, / = l,2,---,iV, (13) 

and 



Replacing the result f lT^ in eq. flT^ . we obtain 



kT A f h 

1 + ^ • (15) 



\J\ 15 V 1^1 

That gives the lowest temperature for which the first excited state is expected to contribute 
to the behavior of the thermodynamic functions. Eq. OlSp shows that T^ax is independent 
of the value of 77,. In this region of parameters of the Hamiltonian ([I]) /(El), the action 
of the crystal field D term is more relevant than the exchange coupling term J and the 
Zeeman term. 

1.2)for -l-A<^<iandA>o U ^ > ^ and A > ^ _ 1 

\^f)B = |l)i®|l)2®---®|0)i8)---®|l)^, / = l,2,---,iV, (16) 

E^-E^ f h D\ 



and 



For this set of parameters t4t and -ryr, eq. (TT2|) gives 



\j\ "^^^ \J\ 

kJ-max 2 /_ h U 



\J\ 15 V \J\ \J\ 



;i8) 



The results (IT^ and fITS]) belong to the same phase B and their difference with respect 
to the first excited state come from the fact that they describe the behavior of the chain 
of spins-1 in distinct ranges of the parameters of the Hamiltonian: in case (1.1) (case 
(1.2)) the effect of Mr is more (less) important than the exchange coupling and the term 
Zeeman effect together. 

2) the first excited state in phase A (see FigJT^): |j| > | and < jjj < jjj — ■^ 

\^?)a = |0)i®|0)2®---®|l)^®---®|0)iv, l = l,2,---,N, (19) 



and 



Ef-E^ f h D , 

2 -— + ^ . (20) 



10 



Replacing the variation of energy ( 120|) in eg. (1121) . we have 

^z™^^if_/L+j^y (21) 

\J\ 15 V 1^1 \J\J 

In the present case the value of T^ax depends on the value of the crystal field D and 
Tmax is a decreasing function of -A. The latter implies that the ground state vector (l8a|) 
is less stable than the vector fl8b|l and fl8c|l under an increasing external magnetic field. In 
phase A of Fig{T^ we have a competition between the crystal field D term and the other 
two terms in the ferromagnetic Hamiltonian ([I])/(l5]). 

In Fig|2]we plot the curves ^^^p ^g ^ function of A for three values of -Mr, covering 

all the regions of -ttt > described in the discussion of the first excited states of the 

ferromagnetic classical spin-1 model (II])/(E])- For tti = 2.5 we have phase A at T = 
and there is a competition between the action of the single-ion anisotropy term and the 
two other terms in the Hamiltonian ([I])/ ([5]). For < -qt < 2 we have a negative slope 
in the curve ^^^p_ This negative slope means that the first excited state can be reached 
without increasing the temperature. 

In FigsJ3^ and 13b we plot the exact curves (solid lines) of A^^ x -p^ with yj-< = 

—2 and 2.5 at rj? = 0.266 and 0.0665, respectively. In FigJS^ the z-component of 
the magnetization of the ferromagnetic model begins at null value with h = and reaches 
its saturated value Aiz = 1 for variations of the magnetic field Ah so that w ~ 7 x 10~^. 

With jjT = 2.5 the function Aiz in FigJSb has a continuous transition from Aiz = to 

Atz = 1 for ^ ~ 4 X 10"'*. The maximum value of the entropy per site for these two 

cases is equal to 4.62 x 10^'^ at h = and -Ar = 2. These are the values of the external 

magnetic field that happens the ferromagnetic phase transitions at T = in diagram [1^. 

In FigjSt we compare the graphs of the z-component of the magnetization of the model 

with Yj\ = 2.5 at two different temperatures: jj-^ = 0.0665 (dotted line) and 0.11 (solid 

line). The latter value of temperature is almost twice the value of ^Tj,"^ with h = 0. It 
is clear from FiglSb and [St the presence of the two plateaus at A^^ = and Aiz = ^ 
at those low temperatures. These plateaus oi Aiz satisfy the Oshikawa, Yamanaka and 
Affleck (OYA) condition[2S]- We remind that this condition, when applied to the periodic 
Hamiltonian ([I])/(I5]) with spin-1, imposes that p{l — Aiz) be an integer, where p is the 
spatial period of the ground state. This is a necessary condition for the occurrence of a 
plateau in the magnetization curve of the one- dimensional spin system[8l [26]. The same 
condition is satisfied in Figj3^. In FigJSb, the transition between the two plateaus happens 
at TjT = 2, the value of the magnetic field per unit of \J\ in which this ferromagnetic 

model suffers a transition between the phases A and B at 77, = 2.5. Note that doubling 

r* I 

the temperature per unit of \J\ does not change much the "step-like" form of the plateaus 
in the curve of A^^. We could say that in the standard ferromagnetic model the plateaus 
are stiff. 

From Figsj3] we verify that outside the transition region tjt where the value of A^^ 
suffers a finite transition, the values of this thermodynamic function correspond to their 
respective value of the ground state in the phase diagram at T = 0. We can use a 
phenomenological approach to fit Aiz in the whole interval of tjt at T < Tmax- We assume 
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that the contributions to Aiz in the transition region of -r^ at T < T^ax come only from 
the ground states in the density matrix operator fITOl) . We have two distinct situations to 
discuss. 

Situation 1): ttt ~ and Mr < ^. In this approximation we obtain for A ^ 0: 



\J\ ~ " "^^^ \J\ — 2- -^^^ "'■'■'■^ c.^i^iw^iiiiu,^x^ii ,,^ ^^^c.xii iwx 1^1 



MMD,h;(3)= _^^, ^ _^^, ^C.c tanh(2/3iV/.). (22) 



Since we consider A to be non-zero and positive we have E^ < Eq . We do not have, a 
priori, a mathematical sound argument to affirm that the approximation (122]) gives the 
right picture of the z-component of the magnetization in the ferromagnetic model at low 
temperatures. 

On the last term on the r.h.s. of result fl22|) we have a product of two limit processes 
in the region of tjt ~ 0: iV — )■ oo (thermodynamic limit) and tjt — >■ 0. The result fl22|) 
depends on the assumption that 

hm AT. A = /(/,; /3), (23) 

'^' 



h 

JJ 



in which the function l{h; (3) is assumed to be finite. The simplest phenomenological 
function l{h; /3) is a linear function 

/(/^;/3) = a-A (24) 

where we take a as a constant. 

Situation 2: nr ~ tji ^ I ^^^ TTf — I- Assuming that for T < T^ax only the ground 
states of the ferromagnetic phases A and B contribute to Aiz , we derive the approximate 
expression of this thermodynamic function Aiz valid for tjt > ■^ at low temperature, 

Mz{J,D,h;{3)^^^^^^^, (25a) 

where 



2|J| 



h f 1 D 



(25b) 



Again in the exponential functions on the r.h.s. of eg. (I25ap we need to calculate the 
product of two limits in the region A ~ ( — ^ + rj])'- iV — )• oo and |4 — )■ ( — ^ + rj])- ^^ 



analogy to eqs. (!23l) and (12^ we write 



lim A^- 



Ar-!-oo 



h / I D 
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Kh;(3). (26) 



In this case we also take the simplest linear function for the phenomenological function 



l(h; (3) = a 



h__( ]_ D_ 

MT^V 2 + M1 



(27) 



where a is a constant. 

In Table [1] we present the values of the parameter a for various values of Mr at distinct 

\^ I 

low temperatures. Some of the values of a were used in the approximate curves of A^^ in 
Figsj3^ and [3]d. 

4 The antiferromagnetic classical spin-1 model 

Fig{T}3 shows the six distinct phases at T = for the AF spin-1 Ising model with a single- 
ion anisotropy term in the presence of a longitudinal external magnetic field. With null 
crystal field (D = 0), the model has only three different phases at T = 0. 

The three terms in the AF (J = -|-1) Hamiltonian ([I])/® compete: the exchange cou- 
pling term J favors neighbor spins to align anti-parallel. The most stable configurations 
of the AF ground state happens when the effect of two terms on each spin in the chain 
are in the same direction. 

We present in appendix [B] the ground state vectors and energies of the six phases of 
the AF Hamiltonian (IT])/([S]) at T = and their respective energy difference between the 
first excited state and its ground state. We also present in that appendix the relation 
between ^'^J"^ and the parameters j and -j such that the condition (ITT|) is valid for the 
classical AF spin-1 model. 

In FigsJH we plot the curves of ^^^^ for various values of y that spans all AF phases 
at T = in the phase diagram FigJT]3. The fundamental difference between Figs. H] and 
m of the ferromagnetic and AF models, respectively, is that the temperature required to 
excite the first excited ferromagnetic state never vanishes. 

In FigJlK with j < we have T^ax = at j = 1. The line 7 = 1 with 7 < separates 
the phases B and G in diagram [Do. The result T^ax = along this line means that the 
classical AF spin-1 model ([I])/® is gap less along the separation of these two phases. 

We also have T^ax = in FigHb with -j = | at j = | and |. These values of j are 
on the lines that separate the phases E ^ G and B ^ E respectively. For the interval 
< 7 < I the AF model ©/dSD is gapless along the lines j = 1 - 7 and j = 1 + 7. At 
J = 2.5 and j = 3.5 with 7 = 2.5 we also have T^ax = in FigJlJo. These values of j 
are on the lines that separate the phases A ^ E and B ^ E respectively. The AF model 
([1])/® is gapless for 7 > ^ along the lines 7 = 7 and j = I + 7. We partially summarize 
the contents of the curves in FigslUby saying that the classical AF spin-1 Hamiltonian 
dl])/® is gapless along the lines that separate the phases in the diagram [T|d at T = 0. 

The presence of plateaus in the thermodynamic functions A^^ of the AF Hamiltonian 
([T]) has been reported previously by some authorsjU [HI HO] for positive values of 7^. In 
the following we discuss the behavior of the thermodynamic functions Ai^ and S of the 
AF version of Hamiltonian ([1]), that is J = 1, by varying the value of the crystal field 
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per unit of J, -j, to span part of the phase diagram of the AF model ([T]) at T = (see 

FigHb). 

Here we draw a more detailed comparison between our phase diagram in FigJTjD with 
the phase diagram of Ref. [8], in which only the case -j > is considered. Before we 
continue the comparison, we should notice that hchen = '^h and Dchen = 2D, where h 
and D are our parameters in the Hamiltonian ([5]). Their Fig. 2b agrees with our diagram 
[Do for -J > 0, except for the fact that in their phase diagram the phase Aiz = m = does 
not distinguish the phases A and G (see the caption of our Fig{T)D) that correspond to 
different ground states. In our phase diagram of the AF model we point out the presence 
of two tri critical points, that is, the points: V{j = 1 and D = 0) and Q{j = ^ and 

For -J < and ^ > we have at T = the phases B and G in diagram [T]d for 
the standard AF model. In this region of parameters of the Hamiltonian we have a 
competition between the states favored by the exchange coupling term and those favored 
by the Zeeman term. 

Figs|5] show the thermodynamic functions Aiz and S versus j with ^ = —2 (long- 
dashed lines). The z-component of the magnetization for jjr < has two plateaus, 
Aiz = and Aiz = 1; at very low temperatures. They satisfy the OYA condition |26]. 
In this region of -j, the transition between these two plateaus happens at j = 1, that 
is the value in which the transition between the phases B and G at T = also occurs, 
for -J < 0. Since along any vertical line in this region of ^ only two phases are crossed 
at T = (see Fig{T]D), the AF spin-1 Ising model ([T]) has only two plateaus at very low 
temperatures for j- < 0. Fig|5K shows the curve Aiz x j at j = 0.266, that is the same 
temperature used in plotting this function for the ferro model in Fig 13^ with -j = —2. 
The curve of Aiz at j = 0.266 still has little resemblance to a "step-like" function. This 
behavior differs from that of the magnetization of the standard ferromagnetic model (see 
FigJS^). By comparing the behavior of the plateaus in the magnetization in the standard 
ferro and the AF models at j = 0.266, we can say that the plateaus of the latter smear 
out even at low temperatures. This happens because although the neighbour spins are 
aligned (anti-aligned) in phase B (G) the exchange coupling term and the Zeeman term 
together with the crystal field D compete in promoting opposing effects on these spins. 
FiglSb presents the entropy per site with -j = —2 at j = 0.266 (long-dashed line). From 
this plot we verify that we cannot approximate the function A^^ by taking into account 
only the contribution of the two ground states of the phases B and G as was done in 
the ferromagnetic version of the model. At j = 0.266, the interval Ah of the magnetic 
field for which the transition between the plateaus Aiz = and Aiz = ^ occurs is such 
that ^ ~ 1. For temperatures of three orders of magnitude lower, the function A^^ as a 
function of j has a "step-like" form and -j^ ~ 10~^. 

Next we consider the AF model ([1]) with jjt = 2.5. In Figj5^ we plot Aiz versus j 

at J = 0.065 (solid line). We verify that at this temperature three plateaus still occur; 
namely, at Atz G {0, |, 1}, for which the OYA condition is satisfied. We have three 
plateaus in this case because the vertical line in the diagram [T]d, localized at 77, = 2.5, 

crosses three phases (B, E and G). The values j = 2.5 and 3.5 where the transitions 
between the plateaus of the z-component of the magnetization occur is the same as that 
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of the transition between the AF phases A ^ E and E ^ B, respectively, of the standard 
model ([T]) at T = 0. Also in this case the curve of Aiz looses its "step-like" form at 
J = 0.266. The width of the transition between the plateaus of Aiz is -^ ~ 0.5 . This 
width reduces to 10^^ for temperatures three orders of magnitude lower than 0.0665 

To understand why the AF function Aiz looses its "step-like" form at j = 0.266 and 
0.665 with -J = —2 and 2.5 respectively, we plot in FigJSjo the function 5 as a function of 
J for these two set of parameters. The function S is non-null around j = 1 with y = —2 
(long-dashed line) and around the values j = 2.5 and 3.5 (solid line). These points in 
the phase diagram [T]d are on the transitions line between the phases A ^ E and E ^ B 
respectively, of the standard model ([I])/([5]) with -j = 2.5 at T = 0. Although the AF 
models are presented at very low temperature, the excited states of the AF model already 
contribute to the thermodynamic function Aiz around the values of j that correspond to 
the lines that separate the AF phases in Fig{T)D. 

The curve Aiz x j with t = | and 2.5 are very similar at low temperatures. For 
y = 0.25, the transition from the plateaus Aiz = to Aiz = \ corresponds, in the 
phase diagram [TJo, to the transition G ^ E, whereas for -j = 2.5 (see FigJS^) the same 
transition corresponds to the phase crossover A :?=i E. The magnetization functions are 
almost insensible if the initial configuration of the spins is either a Neel state or a state 
with all the spins perpendicular to the longitudinal magnetic field. Comparing FigsEti' 
andinb, we see that the specific heat per site, C( J, h, D; /S) [C = — /3^ jda ) distinguishes 

the transitions A ^ E and G ^ E at very low temperatures. In FigJU^ the function C is 
symmetric around -qt = 0.75. 

From Fig. |5^ we could say that the plateaus of A^^ in the AF model ([I])/® are soft, 
in the sense that they are already smeared out at j = 0.26 and 0.0665 with -j = —2 and 
2.5 respectively. 

Since >Vr^( J, h', D'; (3) = Wi(- J, h' , D'; /3), we obtain 

S'"'<^Uh,D;P) = Si-J,h,D;P) (28) 

and 



TWf ^( J, h, D; /3) = Mz{-J, h, D- /3). (29) 

The plateaus in the z-component of the magnetization of the standard ferromagnetic 
(AF) model ([1]) appear in the 2;-component of the staggered magnetization of the stag- 
gered AF (ferromagnetic) model ^. This fact implies that the plateaus in the staggered 
magnetization of Hamiltonian (jl]) satisfy the OYA condition, but for this thermodynamic 
function the AF model has only two plateaus for tjt > 0. 

The behavior of the curves of S^*""^ versus -qt of the staggered ferro and AF models 
is identical to the curves of the entropy per site of the standard AF and ferromagnetic 
models ([I])/® respectively. It is very simple to understand this situation once the phase 
diagram at T = of the staggered ferromagnetic and AF models are depicted by Figs{T]D 
and [1^ respectively. 
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5 Conclusions 

The density matrix method [3] and the /3-expansion of the HFE|2Dj have been apphed to 
the obtainment of the exact expressions of the HFE's of the one- dimensional standard|21] 
and the staggered |22j spin-1 Ising models with {Sz)"^ term in the presence of a longitudinal 
magnetic field. The analytic expressions of these functions have been written in terms 
of exponentials of the parameters of the Hamiltonians (jl]) and (II])/(E])- Our results are 
valid for their respective ferromagnetic (J < 0) and AF (J > 0) models in the interval 
T G (0,oo), extending the results of Ref.[T2] to the AF models. These results do not 
coincide with the results of Ref.[T3]. We present the phase diagram of the standard and 
the staggered spin-models (II])/® and (jl]), respectively, at T = 0. Our result also gives 
the exact HFE of the one-dimensional extended Hubbard model in the ionic limit but in 
the absence of an external magnetic field. 

We have studied the behavior of the z-component of the magnetization Ai^ and the 
entropy S, both per site, of the ferromagnetic and AF models of Hamiltonian (jl])/© at 
low temperatures and their respective staggered versions. 

We have presented the two plateaus of A^^ at low temperatures of the standard fer- 
romagnetic model for tti > I ^^^ show that the value in which the transition between 
plateaus at low temperatures occurs is the same as that of diagram [1^. 

For the AF model ([I])/ ([2]) we show that the number of plateaus of the 2;-component 
of the magnetization at very low temperatures depends on the number of phases of the 
model at T = for a given value of tjt. Our results for the AF model with tji > agree 
with previous results in the literature [TJ [TOl B]- 

The ferromagnetic classical spin-1 model ([I])/® has has highly degenerate excited 
states, but we showed that for temperatures T < T^ax the function Ai^ can be approx- 
imated by the contribution of only two ground states. This fact is not true for the AF 
model because it is gapless for values of j and -j along the line separating the AF phases 
in the diagram [T]d at T = 0. 

The AF specific heat at very low temperatures distinguishes the phase transitions 
E ^ G (at J = 0.75) and E ^ B (at J = 1.25) at j^ = 0.25 and 2.5, respectively, in 
Figsja 

By comparing the plots of the Aiz as a function of j at ^ = 0.266 and 0.0665 to 

y = —2 and 2.5 respectively, we verify that the plateaus of this function of the standard 
ferro and AF models can be called stiff and soft, respectively. 

All the plateaus in the 2;-component of the magnetization of the ferro and AF models 
of the standard and staggered version satisfy the OYA condition. 

As a final comment we should say that the presence of plateaus in the function Aiz 
in the classical spin-1 Ising model with single-ion anisotropy term in the presence of a 
longitudinal magnetic field comes from the stability at low temperature of the ground 
state vector under the action of increasing the norm of the external magnetic field and 
the temperature. 
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FAPEMIG for the partial financial support. M.T.T. also thanks FAPEMIG. M.T.T. 

16 



thanks J. Florencio Jr and J.F. Stilck for enlightening discussions. 

A Calculation of the roots of the third degree 
equation of the spin-1 classical model. 

Krinsky and Furman apphed the transfer matrix method [H [21 [3] to calculate the HFE 
of the classical spin-1 model in Ref.[^. The roots coming from the cubic equations 
have to be real, but they write them in terms of complex quantities. When we handle 
those roots numerically these complex quantities may not disappear. In this appendix we 
recalculate their results for Hamiltonian ([I])/©, showing explicitly that the roots of the 
cubic equation are real. 

Following Ref. p] we obtain that the partition function Zi[J,h,D]P) of the spin-1 
Hamiltonian ([I])/© is equal to 



Z,iJ,h,D;i3) = Tr[V 



Nl 



(A.l) 



where N is the number of sites in the periodic chain and the matrix U for the symmetric 
Hamiltonian i^ is 



V{J,h,D;i3) 



-l3{J+2h+2D) 



-t3{h+D) 



-I3{-.J+2D) 



^-I3{h+D) I ^-i3[-h+D) 

g-/3(-J+2Z)) g-/3(-h+D) g-/3(J-2?i+2D) 



(A.2) 



The matrix \J{J,h,D;P) is hermitian for any value of J, h, D and /3. Its three eigenvalues 
Aj, 2 = 1,2 and 3, are real. The matrices U (see eq. (lA.2|) ) and T (in Ref. [T2]) differ by a 
rearrangement of lines, only. 

In terms of the eigenvalues of U, the partition function (lA.ip becomes: 



Zi(J,/i,D;/3) = (Ai 



,7V 



^< 



N 



Ai 



N 



(A.3) 



in which Ai is assumed to be the eigenvalue of matrix U with the largest modulus. 
In the thermodynamic limit (A^ — )■ oo), the partition function (lA.ip of the model is 



Zi(J,/i,D;/3) = (Ai) 



TV 



(A.4) 



for non-degenerate eigenvalues of U. The HFE of the model is 



Wi(J,/i,D;/3) = -- ln[Ai(J,/i,D;/3)]. 
The eigenvalues Aj, i = 1, 2, 3, are roots of the cubic equation 



(A.5) 



17 



X^ + PX'^ + QX + R = 0, 



(A.6) 



where 



P 

Q 

R 



1 + 2e-''(^+2i?) cosh(2/3/i) = tr[U] 

f3J 



4e 



-2/3Z) 



e '^ cosh(2/3/i) sinh ( ^ ) + 26"^'^^ sinh(2/3J), 



-4/3D 



sinh 



I3J 



sinh(/3J). 



The roots of the cubic equation flA.6p are well known [23| 



(A.7a) 
(A.7b) 

(A.7c) 



where 



with 



and 



Sl 
S3 



-Q cos ( - 



-Q cos 
-Q cos 



P 

+ 3- 
27r\ P 
YJ ^3"' 



27r 

T 



+ 



p 



cos(^) 



R 



-QY 



Q 



3Q + P^ 
9 



(A.8a) 
(A.8b) 
(A.8c) 



(A.9a) 



(A.9b) 



R 



9QP + 27i? + 2P3 

54 



Our previous result Si does not agree with the expression of Xmax of Ref. 
that there are some misprints in their eqs.(18)-(24). 



(A.9c) 
we believe 



B The ground states of the classical AF spin-1 and 



the calculation of T, 



max 



For the AF version of the Hamiltonian ([I])/© we assume that the chain has a even 
number of sites. Let A^ = 2M, in which M is a positive integer. In the thermodynamic 
hmit (A^ — 7- oo) we have M — )■ cxd. 

The ground state vectors at each phase in the diagram [TJd at T = are 

|^o)a = |0)i® |0)2 8)---® |0)2M, (B.la) 

|^o)b = |l)i8)|l)2 8)---®|l)2M, (B.lb) 

|^o)c = |-l)i®|-l)2®---®|-l)y, (B.lc) 

\^q)e = |0)i® 11)2® 10)3® 11)4® ■■■®|0)2A/-i|l)2Af, (B.ld) 

|^o)f = |0)l®|-l)2®|0)3®|-l)4®---®|0)2M-l|-l)2Af, (B.le) 

|*o)g = |l)l® 1-1)2® 11)3® 1-1)4 ®---®|l)2M-l| -1)2M. (B.lf) 

One is reminded that S'f |s)j = s\s)i, s = 0, ±1 and i = 1,2,- ■ ■2M. 

We are interested in the discussion of the thermodynamic function Ai^ at low tem- 
peratures. This function has an even dependence on 4; then we restrict our discussion to 
J > 0. In the following we do not mention the quantum behavior of the phases C and F 
at low temperature. 

The value of the ground state energy of the phases A, B, E and G are, respectively, 



El 

J 
El 

J 

El 
J 

El 
J 



0, 



N 



N 



N 



2D 2h 



D 
7 



■1 + 



h 

1 
2D 
T 



(B.2a) 
(B.2b) 

(B.2c) 

(B.2d) 



We have a much more complex distribution of first excited states Ei in the AF spin-1 
Ising model ([I])/© than in the ferromagnetic version of the model. 

Let us present the difference between the first excited and the ground states of each 
phase A, B, E and G in diagram [TJj. 

1) Phase A: ^ > i and -I < ^. We have 



El - Eq 
J 

which is substituted in eq.( IT2|) to give 



^^ 



J 



15 



h D 

7 + 7 



h D 

j + 7 



(B.3) 



(B.4) 
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2) Phase B: we have two distinct first excited states 

2.1) ^<OandJ>l-^ U f>0andj>l + f. Then 



^^ = 2-1 + ---) (B.5) 



and eq.( 1T^ gives 

k: 

~ J 15 V ' J J 

2.2) ^ < and 1 < ^ < 1 - ^. Then 



kTr^a. 2 / ^^h D ^gg^ 



^^ = 4(-l + ^~), (B.7) 



and eq.( 1T2|) becomes 



^.l(-l4), (B.8) 



3) Phase E: in this phase we have three distinct first excited states. 
3.1) < ^ < i and 1 - ^ < -^ < 1. We have 



with eq.( 1T2|) giving 



Ei-E^ f h D\ ,^ , 



iT™, 2/_^_^fc D. ^g^^j 



r-^ 



J 15 V J J 

3.2) 0<f<iandl<^<l + ^ U 7>iandi + ^<^<l + f. We have 

^^ = 2(1-- + -). (B.ll) 

which gives 

kTmax 2 f h D\ , , 

"""^ - 1 - T + ^ ■ (B-12) 



r-^ 



3.3) |<7<landl<J<i + ^ U 7>landf<$<| + f. We have 
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J \J J 

that substituted in eq.f lT2|) gives 

'^-^ max ^ I t^ ^ 



^.-El _Jh^^ (B.13) 



J 15 V J J 



(B.14) 



4) Phase G: in this phase we have two distinct first excited states. 
4.1) ^ < -1 and < ^ < 1 U -1 < ^ < and ^ < 1 + ^. We have 



Ei-E^ ( h\ 

4 1--. (B.15) 



J \ J 



From eq. flT^ we obtain 



^ . ± (l - ^ I . (B.16) 

4.2) -1 < f < Oand $ < 1 + f U 0<f<iand|<J<l-f U 0<f<iand 
< J < |. We have 



with eq.f ll2p given 



El - E^ ^ 
J ~ 


/ 
1- 

v 


h 
~1~ 


D 

"7 


f^-^max ^ 


1- 


h 


D' 


J '^ 15 


J 


J 



(B.17) 



(B.18) 
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Figure 1: The phase diagrams of the ferromagnetic (J = —1) and the anti-ferromagnetic 
(J = 1) models ([1])/(|2]) at T = 0. (a) The ferromagnetic phases represent the following 
configurations of the 2;-components of the neighbour spins: A — > (0,0), B — > (1, 1) and 
C— ?►(— 1,— 1). (b) The AF phases include the configurations A, B and C that are found 
in the ferromagnetic diagram plus the phases: E — )■ (1, 0), F — )■ (0, —1) and G — )■ (—1, 1). 
This diagram has also two tricritical points: V{A-^ = 1 and D = 0) and Qirj] = 2 ^^^ 
Tjr = |). The phase diagram of the staggered ferromagnetic and the AF models of the 
Hamiltonian (|1]), at T = 0, are given by the figures (b) and (a) respectively. 
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Figure 2: The curves of ,7,''"' for the three distinct first excited states of the ferromagnetic 

model. The dotted line has ^, = —2, the solid line has 171=—^ and the piecewise solid 

line has Mr = 2.5. The vertical dashed line gives the value of the magnetic field where the 
first excited state in the ferromagnetic phase A changes. 
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Figure 3: The exact and approximate expressions of M.^ of the ferromagnetic spin-1 
model ([ID/dSD versus t4. In (a) we have Tjr ^ o „+ A;T 
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2 at f^ = 0.266. The sohd (dotted) 



\J 



hne is its exact (approximate) expression. In (6) we have -fj: = 2.5 at jjr = 0.0665. The 
solid (dotted) line is its exact (approximate) expression. In (c) we compare the exact 
expressions of Aiz with Mr = 2.5 at temperatures pyy = 0.11 (solid line) and 0.0665 



(dotted line). 
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Figure 4: The functions — j^ x j for the AF classical spin-1 model. In (a) we have the 
curves with -j = — 2 (dashed line) and 7 = — | (solid line), and in (&) 7 = | (dashed 
line) and -j = 2.5 (solid line) 
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Figure 5: The thermodynamic functions M.^ and S of the AF spin-1 model ([I])/® are 
plotted versus j with j = — 2 and 2.5. (a) A^. 
IJl = 0.266 (long-dashed line) and with -j 

per site, S, is plotted in (6) for the same set of parameters: -j = — 2 at pjj = 0.266 

(long-dashed line) and -j = 2.5 at |j| = 0.0665 (solid line) 



as a function of j with y = —2 at 
2.5 at ^ = 0.0665 (solid line). The entropy 
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Figure 6: Comparison of the behavior of the function C x j under the transitions G 
and A ^ E at low temperature. In (a) we have the transition G 
(6) is plotted the specific heat per site of the transition A 
curves we have ^ = 0.01. 
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Table 1: The values of the parameter a in eqs. fl21|) and (1271) for several values of 7^7 and 

different values of temperature. -r4 is the interval of variation of the magnetic field in 
which the function M.z varies continuously between A^^ = and A^2 = 1. 
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